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Abstract 

I 
Recently  Van  Loan  and  Demmel  made  conjectures  about  eigenvalue  perturbations.  Van 
Loan's  conjecture  concerned  the  smallest  perturbation  that  makes  a  stable  matrix  unstable, 
and  Demmel's  concerned  the  smallest  perturbations  that  makes  two  matrices  with  disjoint 
spectra  have  a  common  eigenvalue.  We  show  that  the  truth  of  either  of  these  conjectures 
would  imply  the  truth  of  a  third  weaker  conjecture  for  which  we  supply  a  counterexample. 


Recently  Van  Loan  and  Demmel  made  reasonable  sounding  conjectures  about 
eigenvalue  perturbations.  Van  Loan's  conjecture  concerned  the  smallest  perturbation  that 
makes  a  stable  matrix  unstable,  and  Demmel's  concerned  the  smallest  perturbations  that 
makes  two  matrices  with  disjoint  spectra  have  a  common  eigenvalue.  We  show  that  the  truth 
of  either  of  these  conjectures  would  imply  the  truth  of  a  third  weaker  and  equally  reasonable 
sounding  conjecture.  We  then  present  a  counterexample  for  this  third  conjecture  which  hence 
also  serves  as  a  counterexample  to  the  first  two. 

A  stable  matrix  is  a  matrix  all  of  whose  eigenvalues  have  negative  real  parts.  Van  Loan 
[Van  Loan]  recently  made  the  following  conjecture  about  the  smallest  perturbation  of  a  stable 
matrix  A  which  makes  it  unstable  (||-||  denotes  the  2-norm): 

Conjecture  1:  Let  A  be  stable.  Let  B  be  the  closest  unstable  matrix  to  A,  i.e.  B  is  unstable 
and  minimizes  \\A  —  C\\  over  all  unstable  C.  Then  B  has  an  eigenvalue  on  the  imaginary  axis 
with  the  same  imaginary  part  as  some  eigenvalue  of  A. 

If  this  conjecture  were  true,  it  would  lead  to  a  simple  computational  scheme  for 
computing  ||A  -B  ||: 

|[A-B||=    min    cxmin(A  -  I  ImX  /) 

where  o-(A)  is  the  spectrum  of  A  and  i  =  V  —  1. 

sep^(A,fl)  is  the  size  of  the  smallest  perturbations  to  A  and  B  which  makes  them  have  a 
common  eigenvalue: 

sepx(A,fl)  -  min  max(<rmu)(A  -XI)  ,  amin(fl-X/))    . 

Let  a  (A)  denote  the  set  of  eigenvalues  of  A  and  similarly  for  a(B).  Let  co(X)  denote  the 
convex  hull  in  the  complex  plane  of  the  point  set  X.  Demmel  made  the  following  conjecture 
about  the  minimizing  X  in  the  definition  of  sepx: 

Conjecture  2:  The  minimizing  X  in  the  definition  of  sepx  above  lies  in  the  convex  hull 
co({c(A)  ,  a(B)})  of  the  spectra  of  A  and  B. 

If  this  conjecture  were  true,  it  would  greatly  limit  the  region  in  the  X -plane  that  had  to 
be  searched  for  the  minimizing  X. 


. 


.  Bl 

3:"ri1i 


..- 


i 


: 
-6310  nis:>  s#  ,x:  '    ':  '-'J2i<-  ■'- 


zur  . 


I      •  loCJOM   51       (IqOS    UIUOW   I    311 


k« 


■ 


r    I 


+ 


.' 


Of  f:Oj    tl      CObDlbST'lT! 
9lV»  (s.A^   «di  liana 

sis  ?:xii  if  - 

*>   -    Kti  ■       •--      '■ 

I    .     X5V202  ■    ' '    '     ■  - 

•'" 

-,;,--  ai    .'    S 


v.    '.  ■ 


In  this  short  note  we  show  if  either  of  these  two  conjectures  were  true  then  a  third 
weaker  and  equally  reasonable  sounding  conjecture  would  be  true.  Then  we  will  present  a 
counterexample  to  this  third  conjecture.  First  we  need  some  notation.  Define  S(A,e)  as  the 
set  of  all  eigenvalues  of  all  matrices  A  +  8A  for  ||8A  ||^e: 

S(A,t)  =  {X:  det(A  +  8A-X/)  =  0  ,  ||8A||  s  e}    . 

Conjecture  3:  Let  A  be  a  matrix  with  a  single  eigenvalue  X.    Then  S(A  ,c)  is  convex. 

It  is  easy  to  see  how  Conjecture  3  is  implied  by  either  of  the  first  two  conjectures.  First 
consider  Conjecture  1.  Note  that  S(A,e)  is  connected,  since  any  component  must  contain  X. 
As  e  increases,  5(A,c)  grows  from  a  single  point  X  for  e  =  0  to  larger  and  larger  sets.  The 
value  of  €  for  which  this  set  first  touches  the  imaginary  axis  is  the  size  of  the  smallest 
perturbation  that  makes  A  unstable.  Suppose  there  is  a  matrix  A  which  violates  conjecture  3. 
By  multiplying  A  by  a  complex  number  w  of  absolute  value  1  and  adding  a  multiple  a  of  the 
identity,  we  can  rotate  and  shift  the  eigenvalues  of  A  so  that  A  is  stable  and  S(A,e)  makes 
any  angle  to  the  imaginary  axis  we  want.  Since  S(A,(.)  is  nonconvex,  we  can  clearly  choose  co 
and  a  so  that  5(wA+a/,e)  appears  as  in  Figure  1.  By  varying  o>  and  a  slightly  from  these 
values,  we  can  clearly  make  S(a>A  +  a/,e)  either  first  touch  the  imaginary  axis  only  at  a  single 
point  above  the  origin,  or  at  a  single  point  below  the  origin.  Thus  we  can  guarantee  that  it 
does  not  touch  directly  to  the  right  of  the  single  eigenvalue  of  wA  +  al.  Thus  Conjecture  1  is 
clearly  violated.  Therefore  the  truth  of  Conjecture  1  would  imply  the  truth  of  Conjecture  3. 

Now  consider  Conjecture  2.  Again  we  proceed  by  contradiction.  If  Conjecture  3  were 
false,  we  could  find  an  A  with  a  single  eigenvalue  and  an  €  such  that  5(A,e)  were  nonconvex. 
Again  choose  a>  and  a  so  that  S(a>A -fa/.e)  appears  as  in  Figure  1,  with  the  additional 
condition  that  the  two  points  where  S(utA  +  aI,€)  contacts  the  imaginary  axis  are  equidistant 
from  the  origin.  Then  S(-wA  -  o/,«)  is  clearly  the  reflection  of  S(wA  +  c/,t)  in  the 
origin,  as  shown  in  Figure  2.  This  violates  Conjecture  2,  since  the  convex  hull  of  the 
spectrum  of  o>A  +  a/  and  -a>A  -  a/  is  a  line  segment  through  the  origin  passing  between 
the  single  eigenvalue  X  of  wA+a/  and  -X,  and  the  minimizing  X  in  the  definition  of  sepx 
must  lie  at  one  of  the  two  points  of  contact  on  the  imaginary  axis.  Therefore  the  truth  of 
Conjecture  2  would  imply  the  truth  of  Conjecture  3. 

Finally,  we  present  a  counterexample  to  Conjecture  3,  which  is  therefore  also  a 
counterexample  to  Conjectures  1  and  2.    Let 

A  = 

where  fl»l.  A  contour  plot  of  log10(amiD(A  -  X/))  in  the  X  plane  (shapes  of  5(A,e)  for 
various  e)  is  shown  in  Figure  3  (for  B  =  100);  the  nonconvexity  of  the  contours  is  apparent.  In 
fact,  some  of  the  5(A,e)  are  not  even  simply  connected!  From  Figure  3,  we  see  that  0  is  a 
local  maximum  of  the  function  fmjn(A-X).  Thus,  for  example,  S(A,10-3)  is  essentially  a 
disk  with  a  small  hole  near  the  origin.  In  other  words  one  can  make  A  +  8A  have  any 
eigenvalue  in  an  annulus  about  0  with  smaller  ||8A||  than  is  needed  to  make  A  +  8A  have 
eigenvalue  0. 

To  see  how  much  Conjecture  1  can  be  violated,  consider  the  function  tr min(A  —  i\iI), 
where  u,  is  real.  A  plot  of  log10(amin(A  -  /jt/))  versus  ji  is  shown  in  Figure  4  for  B  =  100.  We 
will  show  that  for  \L  =  2~xa  o-min(A -/u./)  is  at  most  33/2  /  (2fl2)  whereas  cmin(A)  is  of  order 
VB ,  which  is  much  larger.  Therefore  S(A,e)  would  touch  the  imaginary  axis  at  about  ±i"2_1/2 
for  e  =  0{B~2)  but  not  the  origin  until  e~fl_1. 

The  proof  is  a  simple  computation. 
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for  fi»l.  When  X  =  i>#0, 

amta(A-/|fc/)  =   ||(A-/m./)-'||-'  : 

(l  +  «p.)3  |n|fl2 

which  as  a  function  of  n,  reaches  its  minimum  33/2  /  (2  fl2)  at  ji  =  2_1/2. 
If  we  let  A  be  n  by  n  and  of  the  same  structure  as  before: 


A  = 


then  o-min(A)  ~  B  '  as  before  and  o-min(A-i>/)  achieves  its  minimum  0(S'  ")  for  \l  =  0(\). 
Thus  for  large  B  and/or  large  n,  using  Conjectures  1  and  2  as  computational  heuristics  can 
lead  to  very  bad  results. 

Note  however  that  the  matrix  A  is  quite  special:  not  only  is  it  defective  but  it  is  nearly 
derogatory.  Of  course  perturbing  A  slightly  would  yield  a  matrix  with  distinct  eigenvalues 
with  similarly  shaped  S(A,(.),  so  defectiveness  per  se  is  not  essential,  but  nearness  to  a 
defective  and  derogatory  matrix.  It  appears  that  if  A  is  far  from  a  derogatory  matrix  (i.e.  A 
can  be  block  diagonalized  with  one  block  per  eigenvalue  using  a  well-conditioned  similarity), 
then  one  cannot  go  too  far  wrong  using  Conjecture  1  as  a  heuristic,  and  since  derogatory 
matrices  are  quite  rare  (in  the  sense  that  a  random  matrix  is  unlikely  to  be  very  close  to  one 
[Demmel]),  the  heuristic  is  likely  to  be  reliable. 
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Figure  3.  Contour  plot  of  log10(o-min(A  -  X/))  in  the  X-plane 
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Figure  4.  Graph  of  log10(amin (A -<>/))  versus  * 
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